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ABSTRACT 


The  design  of  acoustic  transducer  arrays  often  Involves  the 
satisfaction  of  multiple,  conflicting  design  specifications.  The  use 
of  conventional  mathematical  programming  techniques  in  this  endeavor 
have  been  inhibited  by  their  Inability  to  consider  these  multiple 
objectives  within  their  design  formulations.  Most  present  techniques 
formulate  the  design  problem  in  terms  of  a single  objective  subject  to 
a set  of  design  constraints  which  become  difficult,  if  not  Impossible, 
to  solve  if  these  constraints  are  conflicting  in  nature.  Goal 
Programming  is  an  effective  alternative  methodology.  Array  pattern 
synthesis  through  Goal  Programming  provides  the  optimal  solution  to  the 
multiple  criteria  design  of  various  classes  of  acoustic  transducer 
arrays.  The  design  problems  are  formulated  in  terms  of  the  multiple 
objectives  expressed  by  the  system  response  functions  and  the  desired 
characteristics  of  the  design  parameters.  The  solution  of  two 
reasonably  complex  planar  array  designs  is  presented.  The  details  of 
their  formulation  and  solution  provide  an  indication  of  the  ability  of 
this  approach  to  resolve  even  more  complex  array  design  problems. 


CHAPTER  I 


INTRODUCTION 


1 . 1 Purpose  of  the  Study 

This  study  provides,  through  Goal  Programming  [3],  a computerized 
approach  to  the  multiple  criteria  design  of  various  classes  of  acoustic 
transducer  arrays.  Although  the  study  is  concerned  solely  with  the  design 
of  planar  arrays,  the  results  are  Indicative  of  the  ability  of  this 
approach  to  resolve  even  more  complex  array  design  problems. 

The  design  of  acoustic  arrays  is  often  concerned  with  achieving 
a number  of  specific  response  objectives  for  the  main  beam  of  the  array 
while  at  angles  off  this  main  beam  the  pattern  must  lie  below  a specified 
response  level.  Unfortunately,  the  achievement  of  both  of  these  primary 
requirements  is  generally  inhibited  by  their  conflicting  nature.  For 
example,  having  achieved  the  desired  main  beam  response  pattern,  it  is 
relatively  easy  to  decrease  the  response  level  of  the  side  lobes.  This 
is  often  accomplished  by  simply  radiating  a greater  portion  of  the  total 
input  power  from  elements  located  at  the  array  center.  However,  this 
also  produces  a broader  beam,  and  therefore  violates  the  specified 
response  achieved  by  the  main  beam. 

1 . 2 Statement  of  Specific  Objectives 

The  specific  objectives  of  this  thesis  are: 

1.  Develop  a mathematical  model  which  describes  the 
multicriteria  array  design  problem  via  a set  of 
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mathematical  statements  and  extend  this  mode]  to 
the  goal  programming  formulation. 

2.  Achieve  two  specific  acoustic  array  designs  using 
the  Goal  Programming  approach. 

3.  Discuss  the  versatility  of  Coal  Programming  in 
array  design  relative  to  other  tools  used  in  such 
designs. 

4.  Through  the  successful  nature  of  the  results 
achieved  and  the  straightforward  nature  of  the 
problem  formulation,  encourage  the  use  of  Goal 
Programming  in  resolution  of  even  more  complex 
array  design  problems. 

It  should  be  emphasized  here  that  the  goal  of  this  study  is  to 
develop  and  Illustrate  a design  technique,  not  to  provide  specific 
acoustic  array  designs. 


1 . 3 Introductory  Definitions 

Before  analyzing  array  design,  it  is  necessary  to  define  several 
concepts.  The  first  two  definitions  supply  a very  basic  understanding 
of  the  physical  characteristics  involved  in  the  problem  at  hand.  The 
next  three  definitions  provide  an  insight  into  some  of  the  various 
classes  of  acoustic  arrays  in  existence.  Finally,  the  last  two 
definitions  are  provided  to  clarify  the  types  of  mathematical  models 
that  are  employed  to  describe  the  array  design  problem. 

Transducer : 

1.  A transducer  is  a device  for  radiating  or  receiving 
acoustic  waves.  One  of  the  chief  functions  of  the 
transducer  is  to  concentrate  the  radiated  energy 
into  a shaped  beam  which  points  in  the  desired 
direction  in  space  and  to  suppress  the  radiation  in 
other  directions. 

2.  The  radiation  pattern  of  a transducer  is  a graphical 
representation  in  three  dimensions  of  the  radiation 
of  the  transducer  as  a function  of  direction.  The 
radiation  patterns  of  all  practical  transducers 
contain  a main  lobe  and  several  auxiliary  lobes 
called  side  lobes. 
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Arrays : 

3.  A llnejir  transducer  arra^  Is  a set  of  radlatinj; 
elements  arranged  on  a line  as  illustrated  In 
Figure  1.1. 

A.  A planar  transducer  array  Is  a two-dimensional 
distribution  of  elements  located  in  a plane  as 
sI)own  in  Figure  1.2. 

5.  A spatial  transducer  array  is  a three-dimensional 
distribution  of  radiating  elements  as  is  shown 

in  Figure  1.3. 

Mathematical  Models: 

6.  A linear  model  of  a design  problem  describes  all 
objectives  to  be  achieved  by  the  array  elements 
in  terms  of  strictly  linear  functions.  This  is 
often  accomplished  by  transforming  nonlinear 
functions  to  linear  approximations. 

7.  A nonlinear  model  of  a design  problem  describes 
at  least  some  of  the  objectives  to  be  achieved 
by  the  array  elements  in  terms  of  nonlinear 
functions.  Most  array  design  problems  fall  into 
this  category. 

From  these  definitions,  it  should  be  apparent  tliat  there  may  exist 
a nonlinear  model  of  a linear  array,  although  the  expression  itself  seems 
conflicting.  ^ 

1 . 4 The  Chapters  to  Follow 

A number  of  planar  arrays  have  now  been  designed  through  the  Goal 
Programming  approach  with  the  actual  calculations  performed  by  a computer 
algorithm  [3].  The  results  that  have  been  obtained  are  considered 
highly  favorable  and  indicative  of  the  superiority  of  this  approach.  An 
analysis  of  selected,  representative  results  and  the  models  formulated 
to  achieve  them  is  contained  in  the  following  chapters. 
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Chapter  II  presents  a historical  perspective  which  will  briefly 
Introduce  techniques  that  have  been  applied  to  acoustic  array  design. 

This  chapter  Includes  an  Introduction  to  CKial  Prograiranlng  and  a discussion 
on  why  It  was  selected  In  view  of  the  techniques  previously  or  currently 
in  use.  The  formulation  of  the  mathematical  model  of  the  array  design 
problem  and  Its  parallel  Goal  Programming  formulation  are  presented  in 
the  third  chapter.  Finally,  the  results  of  two  array  design  problems 
presented  in  Chapter  IV  provide  an  indication  of  the  power  and  robustness 
of  the  computerized  Goal  Programming  algorithm. 


CHAPTER  II 


HISTORICAL  PERSPECTIVE 

2 . 1 Simulation  Techniques 

In  the  past,  the  chief  means  of  array  design  have  stemmed  from 
studies  using  simulation,  experiments  with  scale  models,  optimization 
methods,  or  some  combination  of  these  three.  The  use  of  simulation  and 
scale  models  often  provides  the  most  costly  and  Inconclusive  results. 

The  results  are  described  as  "Inconclusive"  by  virtue  of  the  fact  that 
these  techniques  can  offer  no  other  criteria  for  the  acceptance  of  their 

I 

results  than  that  the  "optimum"  design  achieved  was  simply  only  the  best 
of  the  relatively  few  examples  tested. 

2 . 2 Classical  Techniques 

As  late  as  1966,  the  best  method  for  the  design  of  linear  arrays 
was  cited  as  the  Dolph-Chebyshev  array  synthesizing  technique  [13].  This 
method  [9]  attempts  to  find  the  array  that  produces  a radiation  pattern 
with  the  narrowest  beamwldth  for  a specified  maximum  side  lobe  level. 
Whereas  synthesis  techniques  employing  Chebyshev  polynomial  expansions 
generally  utilize  a "minimax"  criteria  to  provide  a pattern  whose  deviation 
from  the  desired  pattern  is  minimized,  the  Fourier  integral  approach  to 
-i  synthesis  provides  the  best  fit  in  the  minimum-mean-square-error  sense 

[17].  These  methodologies  Impose  several  constraints  on  the  nature  of 

• the  phasings,  element  spaclngs,  and  excitations  of  the  transducers  which 

« 

(may  Inhibit  tne  achievement  of  a superior  design  [13]. 

I 
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The  Dolph-Chebyshev  method  appeared  during  1946  at  approximate  1 y 
the  same  time  as  the  Woodward-Levins^n  array  synthesizing  technique  [23]. 
The  latter  technique  basically  consists  of  reconstructing  a desired  array 
pattern  from  a finite  number  of  sampled  values.  The  process  provides  a 
pattern  which  achieves  the  desired  pattern  at  a finite  number  of  points, 
but  the  fit  elsewhere  may  be  quite  poor.  The  need  for  a more  systematic 
approach  to  design  problems  gave  rise  to  the  use  of  mathematical 
optimization  techniques. 

2 . 3 Mathematical  Opt Imization  Techniques 

Mathematical  programming  as  it  is  used  in  engineering,  economics 
and  other  fields  involves  the  formulation  of  a mathematical  model  of  the 
problem  and  its  solution  via  one  or  a combination  of  a number  of  optimiza- 
tion techniques.  The  mathematical  model  for  the  design  problem  at  hand 
converts  the  problem  into  a set  of  mathematical  statements  describing 
the  behavior  of  a set  of  response  functions  and  any  restrictions  to  be 
placed  on  the  characteristics  of  the  transducer  array.  Through  the 
response  functions,  the  design  engineer  expresses  his  or  her  objectives 
concerning  the  response  to  be  achieved  within  the  main  beam  and  side 
lobes.  This  limited  discourse  on  mathematical  programming  will  suffice 
as  the  introduction  necessary  to  discuss  the  programming  techniques  that 
have  been  used  to  resolve  array  design  problems.  A more  detailed  analysis 
of  problem  formulation  and  the  process  of  modeling  is  presented  in  the 
following  chapter. 

The  application  of  mathematical  programming  to  Che  design  of 
acoustical  systems  was  little  used  until  the  1960’s.  Since  that  time, 
much  emphasis  has  been  placed  on  the  design  of  linear  models  or  linear 
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approximations  of  the  nonlinear  case.  Although  the  majority  of  array 
design  problems  will  be  nonlinear,  linear  designs  have  been  stressed 
chiefly  because  of  the  Inability  of  most  optimization  techniques  to 
effectively  solve  a design  problem  with  nonlinear  mathematical  statements. 

Linear  programming  techniques  were  employed  by  McMahon,  Hubley, 
and  Moharaned  [15]  in  1972  and  more  recently,  in  early  1976,  by  Wilson  [22] 
to  solve  for  array  designs  which  could  be  expressed  as  a strictly  linear 
set  of  functions.  The  general  linear  programming  problem  involves 
optimization  of  a linear  objective  function  subject  to  both  inequality 
and  equality  constraints  expressed  in  terms  of  the  control  variables. 

The  problem  is  placed  in  a linear  programming  format  and  solved  by  the 
simplex  method  [2].  The  approach  of  McMahon  et.  al.  is  typical  of  the 
linear  programming  solution  techniques.  Within  the  linear  programming 
format,  they  establish  a minimizing  objective  function  which,  through  its 
cost  coefficients,  ranks  the  satisfaction  of  the  design  constraints  in 
this  order:  (1)  peak  response  constraints,  (2)  main  beam  offpeak 

response  constraints  and  (3)  constraints  defining  the  desired  response 
in  the  minor  lobe  region.  Lasdon  et.  al.  take  the  linear  array  problem 
one  step  further  and  apply  nonlinear  techniques  to  the  design  of  linear 
arrays  [13,21]. 

Acknowledging  the  nonlinear  nature  of  most  array  design  problems, 
Lasdon,  Waren  and  Suchman  focus  on  the  study  of  nonlinear  acoustic 
transducer  arrays  [14]  and  solve  the  nonlinear,  constrained  optimization 
problem  using  the  Sequential  Unconstrained  Minimization  Technique  (SUMT) 
of  Flacco  and  McCormick  [1].  The  problem  formulation  in  this  case  consists 
of  the  minimization  of  a linear  function  which  is  an  indication  of  the 


amount  by  which  the  solution  vector  fails  to  meet  the  design  specifica- 
tions supplied  by  the  nonlinear  constraints.  The  transducer  design 
formulation  then  involves  the  minimization  of  a function  which  reflects 
the  underachievement  of  the  array  design  specifications  supplied  by  a 
set  of  nonlinear  constraints  which  define  the  desired  response  levels  of 
the  main  beam  and  side  lobes.  Responding  to  the  difficulty  inherent  in 
the  solution  of  optimization  problems  involving  nonlinear  constraints, 
the  SUMT  technique  uses  a penalty  function  to  convert  an  Inequality 
constrained  problem  to  a sequence  of  unconstrained  minimization  problems. 
If  the  first  derivatives  of  this  sequence  of  functions  exist  and  can  be 
easily  computed  from  analytic  expressions,  then  Lasdon  et.  al.  use  the 
descent  method  of  Fletcher  and  Powell  [10]  for  minimization.  If  not, 
then  modified  methods  developed  by  Powell  [16]  or  Stewart  [18]  (which 
solve  for  the  minimum  without  calculating  derivatives)  are  employed. 

The  SUMT  technique  in  conjunction  with  the  unrestricted  methods  of 
Powell  and  Stewart  is  not  a particularly  powerful  one  as  it  provides 
only  the  approximate  location  of  the  optimum,  which  may  not  be  very 
precise  when  this  lies  In  a sharp  comer  [1].  (Sharp  ridges  are 
typically  found  on  the  transformed  surfaces  created  by  the  Interior 
Penalty  Method  used.)  Unrestricted  optimization  techniques  such  as 
Powell's  are  also  known  to  perform  poorly  computationally  with  a large 
number  of  variables  [11].  In  addition,  the  SUMT  technique  Is  not  a very 
robust  one,  requiring  a feasible  starting  point  from  which  to  begin  the 
search  for  an  optimum.  Although  there  exists  techniques  such  as  those 
of  Box,  Carroll  or  Hensgen  [1]  to  find  a feasible  base  point,  these  often 
involve  a computationally  lengthy.  Iterative  procedure. 
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Given  the  multimodal  nature  of  most  design  problems.  Inherent  in 
all  nonlinear  problems  Is  the  lack  of  assurance  that  a local  optimum  Is 
the  global  value.  Generally,  It  Is  necessary  to  start  from  different 
points  in  pursuit  of  the  global  value  and  In  order  to  avoid  acceptance 
of  a poor  local  optimum.  However,  with  the  extensive  set  of  restrictions 
generally  present  in  array  design  problems,  it  is  difficult  to  obtain 
valid  alternative  starting  points  which  satisfy  all  the  constraints  and 
are  significantly  different. 

Avoiding  the  need  for  an  initial  feasible  solution,  H.  S.  C.  Wang 
employed  the  method  of  Lagrange  multiplier  with  main  beam  normalization 
as  a single  constraint  in  an  attempt  to  minimize  the  side  lobe  power  [19]. 
Assuming  that  the  output  power  is  a quadratic  of  the  amplitude  shading 
coefficients,  this  method  leads  to  a set  of  linear  simultaneous  equations 
which  he  solves  using  a standard  numerical  technique.  Essential  in  the 
use  of  this  method  is  the  existence  of  gradients  that  can  be  found 
directly  from  analytical  expressions. 

The  method  of  Lagranglan  multipliers  is  not  commonly  used  in 
optimizing  very  complex  systems  as  the  required  calculations  become 
intolerable.  Further,  the  final  product  of  such  calculations  is  only  a 
stationary  point  and  it  must  be  evaluated  to  guarantee  that  it  is  an 
optimum  [1]. 

In  more  recent  work,  Wang  uses  the  method  of  steepest  ascent  in 
pattern  synthesis  [20].  His  problem  formulation  consists  of  objective 
functions  which  are  "suitable"  linear  combinations  of  the  acoustical 
power  radiated.  The  development  of  the  mathematical  mtidel  of  the  problem 
is  an  iterative  procedure  with  the  "suitable"  linear  combinations  of  the 
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response  functions  In  the  objectives  derived  from  a number  of  attempts 
at  the  "appropriate"  combinations.  Solving  for  these  objective  functions, 
Wang  employs  the  method  of  steepest  ascent  to  search  for  the  character- 
istics of  the  array  which  give  a maximum  value  of  an  objective  function. 

The  method  of  steepest  ascent  is  basically  a cyclic  search 
technique  which  approximates  the  line  of  steepest  ascent.  At  each  point 
in  the  search,  the  local  function  gradient  is  evaluated  and  search  for 
the  restricted  optimum  is  made  along  the  resulting  direction  of  steepest 
ascent.  Although  each  one-dimensional  search  is  searching  in  the  "best" 
direction,  this  method  is  not  very  effective.  It  can  easily  be  seen  that 
the  direction  of  steepest  ascent  as  it  is  defined  here  is  only  a local 
property  which  does  not  account  for  interactions  between  variables  [Ij. 
Computationally,  the  technique  is  cumbersome  since  the  derivatives  of  the 
objective  function  must  be  evaluated  at  each  base  point  in  the  search. 

As  in  the  method  of  Lagrange  multipliers,  the  existence  of  easily 
established  gradients  is  critical  in  the  use  of  this  technique.  If  the 
gradients  cannot  be  found  directly,  they  must  be  estimated  by  local 
explorations  which  may  be  quite  excessive. 

2 . 4 Goal  Programming 

Since  unrestricted  methods  of  optimization  are  generally  better 
than  the  constrained  optimization  techniques  [1],  techniques  such  as  SUMT 
have  focused  on  the  replacement  of  a constrained  optimization  formulation 
by  an  equivalent  problem  without  constraints.  Goal  Programming  [3],  a 
tool  for  multicriteria  design,  was  developed  by  Charnes  and  Cooper  for 
the  strictly  linear  problem  in  the  1950's  [8]  and  has  since  been  expanded 
by  Ignlzio  [3]  to  consider  the  nonlinear  model  as  well.  The  Goal 


Programming  problem  formulation  expresses  the  array  design  problem  In 


terms  of  its  multiple,  conflicting  objectives  (such  as  the  desired  main 
beam  and  side  lobe  response  levels).  The  general  formulation  establishes 
the  problem  as  one  of  minimizing  the  difference  between  the  pattern 
achieved  by  the  antenna  array  design  and  the  desired  response  levels 
expressed  in  the  objectives.  The  specific  Goal  Programming  method  used 
in  this  thesis  to  achieve  the  optimum  array  design  is  a modification,  by 
Ignizio,  of  the  numerical  search  technique  of  Hooke  and  Jeeves  [12].  This 
method,  also  known  as  modified  pattern  search,  begins  with  an  initial 
base  point  and  perturbs  about  this  point  along  all  coordinates  of  the 
array  in  search  of  an  improved  solution  in  terms  of  the  prioritized 
design  objectives.  The  Hooke  and  Jeeves  method  is  referred  to  as  an 
accelerated  search  technique  [1]  as  it  evaluates  the  success  along  any 
direction  and  accelerates  the  ongoing  search  along  successful  paths  and 
reduces  the  distance  explored  in  the  event  of  failure  to  improve  the 
design.  Further  discussion  of  the  modified  pattern  search  algorithm  can 
be  found  in  Appendix  B. 

The  modified  Hooke  and  Jeeves  optimization  technique  was  employed 
in  lieu  of  a number  of  numerical  search  techniques  as  it  is  a simple  and 
effective  method  of  minimization  which  is  able  to  handle  a large  nuraoer 
of  variables  and  multiple  objectives  and  possesses  a range  of  convergence 
criteria.  Numerical  search  methods  require  only  the  evaluation  of  the 
achievement  function  at  a particular  location,  whereas  gradient  techniques 
require  both  function  and  gradient  values  to  be  found  at  each  location; 
if  the  gradient  cannot  be  found  directly  from  an  analytical  expression, 
it  must  be  estimated  by  local  exploration  demanding  several  experiments. 


Therefore,  in  light  of  the  size  of  most  design  problems,  pattern  search 
is  much  more  efficient  in  terms  of  computation  time. 

Optimization  using  nonlinear  Goal  Programming  does  not  require  a 
feasible  starting  point,  one  satisfying  all  objectives,  as  do  techniques 
such  as  SUMT  and  Complex  [1],  removing  tlie  initial  task  of  seeking  what 
may  be  a nonexistent  point.  Therefore,  this  method  is  more  efficient  in 
its  search  for  the  global  optimum  through  its  ease  in  choosing  significantly 
different  starting  points.  However,  as  in  all  methods  of  nonlinear 
optimization,  there  Is  no  known  way  to  guarantee  that  the  solution  or 
system  design  is  the  global  optimum.  For  these  reasons,  the  modified 
pattern  search  method  of  nonlinear  Goal  Programming  was  selected  as  the 
tool  for  acoustic  array  design  In  this  study. 


CHAPTER  III 


PROBLEM  FORMULATION 

3.1  Mathematical  Model  Formulation  (Planar  Arra^^ 

The  Initial  step  In  the  Goal  Programming  formulation  of  array 
design  Is  the  development  of  a preliminary  matliemat  ica  I model  wliich 
essentially  describes  the  multlcrlterla  design  problem  via  a set  of 
mathematical  statements.  This  Introduction  to  mathematical  models  of 
planar  arrays  will  be  followed  by  the  extension  of  the  preliminary 
mathematical  model  to  the  Goal  Programming  formulation.  A specific 
planar  array  design  model  formulation  will  be  used  to  Illustrate  the 
process . 

All  planar  array  design  problems  studied  were  symmetric  about 
the  horizontal  and  vertical  axes  In  the  plane  and  uniformly  spaced 
at  one-half  wavelength  apart.  However,  after  reviewing  the  Goal 
Programming  formulation  which  follows.  It  should  be  apparent  that  this 
approach  can  easily  be  extended  to  more  complex  array  design  problems. 

A primary  goal  of  planar  array  design  Is  to  form  a certain 
radiation  pattern.  The  associated  design  objectives  are  expressed  In 
terms  of  a response  function  or  set  of  response  functions  which  describe 
the  gain  or  directivity  of  the  array.  This  gain  Is  measured  over  the 
standard  6,4)  coordinate  system  (In  which  0 Is  the  angle  of  incidence 
as  measured  from  the  positive  Z axis  and  4>  is  the  roll  angle  from  the 
X axis  of  the  projection  In  the  X-Y  plane,  (see  Figure  3.1)  and  will  be 
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denoted  as  g(0,4')  • Kach  specific  value  of  (G,(}))  considered  in  the 

objectives  is  thus  assigned  a design  response,  g^(0,(}))  . 

The  response  objectives  of  the  particular  design  problems 

studied  were  to  obtain  a main  beam  with  uniform  response  and  to  have  low 

response  elsewhere.  Defining  index  sets  I and  I such  that  I 

A B A 

contains  those  values  of  1 for  which  g^(0,<t))  must  be  uniform  at  a 
specified  level  Y and  I those  for  which  g.(0,ji)  is  to  lie  below 

D 1 

another  level  B , then  the  design  objectives  may  be  expressed  as: 


g,  (9.4>)  = Y 1 e I. 

1 A 


(3.1) 


and 


gj^(0,4>)  16  i e Ig  . 


(3.2) 


Elaborating  on  the  response  functions,  the  power  response  objectives  of 
a symmetric  planar  array  follow  (see  Appendix  A for  the  derivation  of 
these  functions  from  the  general  case): 


gj(0.<l'^ 


= 10  log 


M 


M 2 

16[(  E A • cos  a • cos  B • cos  C ) + 

, m m m m 

m=l 


( E A • sin  a cos  B • cos  C ) 
.,  m m w TO 

m=l 


= Y 


where  lei. 


(3.3) 


gj(0.4') 


= 10  log 
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M 2 

16[(  E A ■ cos  a • cos  B • cos  C ) + 

, m m m m 
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( E A • sin  a • cos  B • cos  C ) J < B 
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where  i c I 
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(3.4) 


given  that  : 
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cos  <p) 


(3.5) 


and 
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m tn 


and 


-j-  (Y  • sin  0 • sin  4>) 
A m 


(3.6) 


the  desired  response  level  for  the  main  beam, 

the  desired  response  level  for  the  side  lobes, 

the  set  of  responses  pertaining  to  the  main  beam, 

the  set  of  responses  pertaining  to  the  side  lobes, 

the  number  of  sets  (of  four)  array  elements  in 
quadrantal  symmetry, 

the  operating  wavelength  of  the  array, 

the  cartesian  coordinates  (in  wavelengths)  of  the  m*"^ 
set  of  four  elements, 

the  phase  angle  associated  with  the  m*’^  set  of  four 
elements , 

the  shading  factor  of  the  m^^  set  of  elements. 


The  particular  design  problems  studied  Involved  obtaining  a main 
beam  with  uniform  response,  within  3 dB  if  possible,  and  a low  response 
elsewhere  (a  goal  of  -20  dB  was  set).  The  control  variables  which  were 
permitted  to  be  regulated  to  achieve  the  specified  response  levels,  were 
either  the  phasing  of  the  array  elements,  their  amplitude  or  both. 
Therefore,  the  design  problem  becomes  that  of  finding  the  vector  a.  A, 
or  (a,  A)  best  sat  isfylng  Equations  (3.3)  and  (3.4).  Due  to  the 
symmetry  Inherent  in  the  problem,  it  is  only  necessary  to  specify  the 
elements  of  a single  quadrant.  The  design  objectives  can  now  be  restated 


as : 
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g^(0,<t>)  - r(0“,0“)  = 0 dB  1 C 

and 

g^(0,4))  - g(0“,0°)  < -20  dB  1 E 

Each  design  objective  of  Equations  (3.7)  and  (3.8)  corresponds  to 
the  desired  radiation  level  at  some  value  of  0 and  (J)  . The  size  of 
and  Ig  , the  number  of  objectives  constraining  the  main  beam  and 
side  lobe  levels  respectively.  Is  determined  by  the  number  of  values  of 
9 and  ({)  that  must  be  considered  so  as  to  attain  the  desired  pattern. 
Oftentimes,  If  the  range  between  values  of  6 and  <})  defined  In  the 
model's  objectives  Is  "too  large"  (Implying  that  a small  number  of 
radiation  levels  have  been  Included  In  the  model) , the  pattern  will  attain 
the  desired  level  within  the  constrained  area,  but  with  violations  of  the 
desired  response  occurring  within  those  regions  not  specified  In  the 
objectives.  On  the  other  hand,  it  Is  essential  to  keep  the  number  of 
objectives  at  a minimum  as  computation  time  and  computer  storage 
requirements  rapidly  Increase  as  the  number  of  objectives  Increases. 

Due  to  the  symmetry  of  the  models  considered,  the  range  on  0 and 
(})  Included  In  the  response  objectives  can  be  reduced  to  extremes  of  0° 
and  90°  (l.e.,  the  first  quadrant).  It  Is  Important  to  Incorporate 

Into  the  model  any  properties  of  the  design  that  can  reduce  the  Inherent 
~ size  of  the  problem. 

3 . 2 Nonlinear  Goal  Programming  Model  Formulation 

The  second  phase  of  the  problem  formulation  Involves  the  extension 
of  the  preliminary  mathematical  model  to  the  goal  programming  formulation. 
This  Involves  the  addition  of  deviation  variables  In  the  objective 

i 

f 

k 


(3.7) 

(3.8) 
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functions  and  the  statement  of  the  "achievement  function"  In  terms  of 
these  variables. 

On  the  left-hand  side  of  the  goal  programming  formulation,  each 
objective  function  will  contain  a negative  and  positive  deviation  variable 
(expressed  as  n^  and  p^  respectively).  The  value  of  n^  reflects 
the  negative  deviation  from  the  desired  response  level  for  any  specific 
solution  vector,  while  the  value  of  p^  reflects  the  positive  deviation. 
Letting  b^  equal  the  desired  response  level  for  a particular  g^(6,(p)  , 
every  response  objective  function  can  be  expressed  as  an  equality  of  the 
form: 


1 

I I 


Looking  at  the  side  lobe  objectives  for  clarification  of  the  role 
of  the  deviation  variables,  if  the  radiation  level  of  a particular  gj^(6,<f>) 
is  below  that  required  20  dB  dropoff  from  the  main  beam  level,  the 
negative  deviation  variable  (n^)  will  take  on  the  value  of  the  additional 
difference.  Conversely,  if  the  value  of  the  response  function  Is  greater 
than  the  design  specifications,  the  positive  deviation  variable  (p^^)  will 
assume  the  value  of  the  design  violation. 

An  expression  of  the  design  objectives  in  terras  of  the  deviation 
variables  follows: 


- g(0“,0')  + n^  - p^  = 0 dB  i e (3.9) 

and 

gj(0,4>)  - g(0“,  0")  + n^  - p^  » -20  dB  1 G . (3.10) 
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A key  assumption  of  nonlinear  Goal  Programming  is  that  the  design 
engineer  can  establish  preemptive  priorities  [4]  for  each  objective  or 
group  of  objectives.  The  first  priority  is  always  assigned  to  any 
absolute  objectlve(s) , l.e.,  any  objectlve(s)  whose  satisfaction  is 
absolutely  necessary  for  the  resulting  solution  to  be  Implementab le . 

Next,  priorities  are  associated  with  all  remaining,  nonabsolute  objectives. 

In  solution,  nonlinear  Goal  Programming  initially  satisfies,  as 
best  as  is  possible,  the  objectives  with  the  highest  priority.  Continuing 
on  to  the  next  priority  level,  the  objectives  at  that  level  are  satisfied 
as  closely  as  is  permissible  without  degrading  the  achievement  of  any 
higher  priority  objectives.  This  process  continues  until  all  priority 
levels  have  been  considered  and  resolved  as  well  as  possible.  It  can  be 
seen  from  this  process  that  objectives  at  a higher  priority  level  are 
immeasurably  more  important  than  those  of  a lower  level. 

Judgmental  or  inherent  weights  may  also  be  assigned  to  the 
commeasurable  objectives  or  goals  within  a priority  level.  These  weights 
must  be  nonnegative  numbers.  When  trying  to  achieve  a difficult  response 
pattern  in  antenna  design,  the  weighting  factors  may  be  assigned  in  such 
a way  that  the  algorithm  will  strive  to  achieve  the  desired  pattern 
within  an  area  of  resistance  at  the  expense  of  the  remaining  response. 

The  determination  of  the  "correct"  weighting  factors  for  a given  design 
is  often  an  iterative  process  with  completion  of  the  process  only  upon 
arrival  at  a "satisfactory"  antenna  array  design. 

In  order  to  remove  any  ambiguity  concerning  the  deviation  variables 
in  the  forthcoming  achievement  function,  it  is  necessary  to  present 
another  restatement  of  the  problem  formulation: 
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gj(0,4))  - g(0“,0“)  + = 0 dB  1 ^ 

and 

g,(0.(J>)  - g(0°,0°)  + n - p = -20  dB  j t U • (3.12) 

J J J B 


Therefore,  the  single  goal  of  the  design  problem  as  defined, 
without  any  restrictions  on  the  range  of  the  phasing  or  shading  coefficients. 
Is  to  minimize  the  quantity  (n^  + p^) , which  demands  a uniform  response 
in  the  main  beam,  and  p^  , which  requires  a response  In  the  side  lobes 
less  than  or  equal  to  a level  20  dB  below  that  of  the  main  beam.  This 
goal  is  expressed: 


gj(n,p)  = (n^  + p^  + Pj)  i C 1^  , j C Ig  • 

The  formulation  of  the  achievement  function  is  the  final  step  In 
the  statement  of  the  Goal  Programming  model.  The  achievement  function  Is 
a linear  function  of  the  deviation  variables  which  will  always  be 
minimized.  In  order  to  maintain  the  preemptive  priorities  of  each  goal, 
the  achievement  function  Is  always  an  ordered  vector  of  the  form  (3): 


where 


and 


a = (gj^(n,p),  g^Cn  n) , . . gj.(n,p)}  , 

(1)  gj^(n,p)  Is  a linear  function  of  the  deviation  variables, 

(2)  k Is  the  priority  l^evel  associated  with  the  objectives 
expressed  In  g|^(n,p), 

(3)  the  dimension  of  a represents  the  number  (K)  of 
preemptive  prlorit.y  levels  among  the  objectives, 

(A)  the  number  of  preemptive  priorities  are  equal  to  or  less 
than  the  total  number  of  objectives, 

(5)  a will  always  be  minimized. 
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The  achievement  function  for  the  design  problem  formulated  to 
this  point  Is: 


minimize 


a = {gj^(n,  p)} 


or 


minimize  a = 


= Un^  Pi  • 


This  achievement  function  assumes  equally  weighted  priorities.  If  a 
weight  of  w^  were  associated  with  each  of  the  objectives  of  the  main 
beam,  gj^(n,  p)  would  read  [w^(n^  + p^)  Pj ) • 

In  final  solution,  the  achievement  vector  will  be  equivalent  to 
0 if  all  of  the  objectives  were  met  and  will  take  on  positive  values  if 
all  could  not  be  met.  The  solution  of  a design  model  is  considered 
optimal  if  the  corresponding  value  of  a is  the  same  or  preferred  to 
the  value  of  a for  any  other  feasible  solution. 

The  final,  completed  formulation  of  the  planar  array  design  model 
follows : 

Find  the  solution  vector  (phasings,  amplitudes  or  both)  so  as  to 
minimize : 


I 


a = {(n^  Pi  Pj^^  ’ 

such  that : 


gj^(0.4>) 

- g(0“,0“)  + nj  - “ 0 0B 

' ^ 'a' 

(3.13) 

gj  (0.4)) 

- g(0“,0“)  + - Pj  “ -20  0B 

^ 'b' 

(3.14) 

and 

n,  p ^ 0 
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We  now  introduce  the  formulation  of  a specific  antenna  design 
problem.  The  problem  to  be  analyzed  is  the  design  of  a 12  x 12  circular, 
symmetric,  uniformly  spaced  planar  array.  The  square  grid  has  the  corners 
truncated  so  that  the  array  approximates  a circular  aperture.  In  this 
analysis  both  the  amplitude  and  phasing  will  be  variable.  Therefore, 
the  example  to  be  analyzed  seeks  the  optimum  phasing  and  shading 
coefficients  of  the  104  element  array.  Due  to  the  symmetry  iniierent  in 
the  problem,  it  is  necessary  to  consider  only  four  sets  of  26  elements. 
(Figure  3.2  indicates  the  location  of  the  26  elements  in  the  first  quadrant 
of  the  X-Y  plane.) 

The  pattern  to  be  achieved  is  a circular  30°  diameter  beam.  Thus, 
the  design  specifications  require  a circular  30°  beam  with  uniform 
response  and  all  angles  off  the  main  beam  at  a level  20  dB  below  the  main 
beam  response.  In  addition  to  these  (ideal)  design  objectives,  it  is 
also  desirable  to  constrain  the  characteristics  of  the  array  elements  to 
allow  for  easy  implementation.  The  phase  angles  are  defined  within  the 
range  of  0°  to  360°  and  the  shading  coefficients  are  permitted  to  take 
on  only  positive  values. 

The  formulation  of  the  goal  programming  model  established  to 
satisfy  these  design  specifications  follows: 

Find  the  vector  (Aj^,A2»  • • •,  A^^,  , a^.  ....  ot^^)  so  as  to 

minimize : 

a = {(nj^  + n^  + . • • + + \>2i  + • . • + P52  ”53 

. . . + ("3  ^1  ’ 


.*‘22  hi  hi 


Figure  3 


.2.  Array  Elements  In  the  First  Quadrant  of  the  X-Y  Plane 
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such  that: 


a + n,  - 
m k 

Pk 

= 0° 

k = 

1,2,  . . 

.,26, 

(3.16) 

a + n,  - 
m 1 

Pi 

= 360° 

1 = 

27,  . . 

..  52, 

(3.17) 

A + n 
m n 

Pn 

= 0 

n = 

53,  . . 

..  72, 

(3.18) 

m = 

1.2,  . . 

.,  26 

- 

g(0° 

.0°)  + - 

Pf  = 

0 dB 

1 e 1 ' 
A 

(3.19) 

theta  ranging 

from 

0°  to  15 

° for 

values  of 

phi  from  0° 

to  90° 

gj(6,4>)  - 

g(0° 

,0°)  + - 

-20  dB 

^ ^ ^b' 

(3.20) 

for  theta  ranging  from  16°  to  90°  at  all  values  of  phi  from  0°  to 
90° , where : 


g^(0.4>)  = 10  log 


r 

16((  E 

■-  m=l 


A • cos  a • cos  B • cos  C ) + 

m m m m 


( E A • sin  a • cos  B • cos  C ) 
, m m m m 

m=l 


and  n,  p ^ 0 

Aga in , 


B = (x  • sin  0 ■ cos  4>) 
m A in 


C * (y  • sin  6 • sin  (|>) 

in  Am 


2b 

Note  that : 

1.  The  size  of  I ' and  I ' Is  dependent  upon  the  size  of 
the  sampling  of  radiation  levels  necessary  to  attain  the  desired  pattern 
without  response  violations  within  the  undefined  regions. 

2.  The  achievement  vector  as  it  is  written  has  equally 
weighted  objectives  which  assumes  the  design  achievement  will  not  be 
confronted  with  resistance  in  any  particular  area.  This  will  seldom, 
if  ever,  be  the  case  in  practice. 

3.  The  operating  wavelength,  A , for  all  of  the  design  prob- 
lems is  assumed  to  be  1. 

4.  The  cartesian  coordinates  (x  , y ) are  expressed  in 

m m 

intervals  of  1/2  wavelength. 

The  subroutine  (YVALUE)  required  to  be  supplied  by  the  user  of 
the  general  computerized  Goal  Programming  code  to  define  this  particular 
model  formulation  for  the  optimization  algorithm  is  presented  in  Appendix 
C.  In  addition.  Appendix  B provides  further  explanation  on  Goal 
Programming  and  the  design  problem  formulation. 

In  the  next  chapter,  we  present  the  results  of  the  application 
of  the  nonlinear  Goal  Programming  code  to  both  this  problem  (l.e.,  a 
circular  beam)  and  also  to  the  problem  of  forming  a rectangular  beam. 
Although  numerous  other  problems  were  solved,  these  two  problems  are 
typical  of  the  approach  and  results  obtained. 


CHAPTER  IV 


EMI’IRICAL  RESULTS  | 

4 . 1 Int  roduct ion  1 

The  solutions  to  all  array  designs  discussed  herein  were  obtained 
using  the  nonlinear  Goal  Programming  code  as  listed  and  documented  In  the 
appendix  of  Reference  3.  All  computer  analyses  were  performed  by  an  IBM 
370. 

The  12  X 12  array  designs  presented  herein  are  the  result  of  a 
three  stage  experiment.  Initially,  the  design  problems  were  formulated 
to  solve  for  a number  of  6 x 6 planar  arrays.  The  study  was  then 
extended  to  8 x 8 planar  arrays.  Finally,  the  12  x 12  models  were 
resolved . 

The  study  began  with  the  smaller  arrays  for  several  reasons. 

First,  the  relatively  small  solution  times  permitted  extensive  manipula- 
tion of  the  parameters  of  the  goal  programming  formulation  (weighting 
factors,  number  of  response  objectives,  etc.)  at  a far  smaller  cost  (the 
computation  time  at  each  iteration  of  the  12  x 12  array  problem  was 
approxlm  "ely  five  times  greater  than  its  6x6  counterpart).  These 
manipulations  provided  an  indication  of  the  responsiveness  of  the  design 
models  to  changes  in  the  parametric  factors.  Second,  the  analysis  of 
the  smaller  models  provided  trade-off  policy  concerning  the  computation 
time  involved,  the  size  of  the  problems  and  the  number  of  Iterations 
till  solution.  Third,  any  direct  comparison  of  the  Goal  Programming 


results  with  those  of  more  conventional  methods  has  to  be  performed  on 
a small  scale  as  most  conventional  methods  have  been  unable  to  consider 
the  larger  problems  that  have  been  resolved  through  Goal  Programming. 

The  results  of  this  three  stage  experiment  in  planar  array  design 
may  be  summarized  as  follows: 

1.  The  patterns  achieved  compare  quite  favorably  with 
the  results  obtained  by  other  conventional  approaches. 

This  Is  an  Impressive  achievement  In  view  of  the 

much  more  extensive  and  costly  nature  of  the  research 
done  using  conventional  methods. 

2.  When  compared  with  more  traditional  optimization 
methods,  the  Goal  Programming  method  has  achieved 
much  better  solution  times  in  relatively  small 
problems.  If  the  traditional  methods  were  extended 
to  handle  larger,  more  complex  problems,  this 
advantage  should  escalate  in  favor  of  Goal  Programming. 

3.  Although  Goal  Programming  can  solve  large  problems 
(12  X 12  array  formulations  involved  52  variables 
and  up  to  193  objectives),  it  is  suggested  that 
those  using  the  technique  begin  with  a relatively 
simple  problem  in  order  to  become  acquainted  with 
the  technique  and  its  options.  This  reduces  the 
computation  time  and  costs  in  the  developmental 
stages  of  the  experiment. 
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A. 2 Design  Specifications 

Array  patterns  are  often  classified  according  to  their  shape.  Two 
shapes  often  sought  in  practice  are  circular  beam  patterns  and  rectangular 
patterns.  The  circular  beam  attempted  here  was  to  have  a 30°  diameter, 
i.e.,  the  3 dB  beamwldth  about  the  beam  axis.  The  rectangular  beam,  a 
more  difficult  pattern  to  realize  in  practice,  attempted  to  achieve 
72°  X 36°  dimensions.  Both  of  these  design  problems  were  to  obtain  the 
specified  circular  or  rectangular  beam  with  uniform  main  beam  response 
(within  3 dB  if  possible)  and  low  response  elsewhere  (20  dB  below  that  of 
the  main  beam).  The  model  formulation  for  the  circular  pattern  as  such 
was  presented  in  Chapter  III. 

A . 3 Beam  A Circular  Array  Pattern 

The  amplitude  and  phasing  ascertained  through  the  nonlinear  Goal 
Programming  approach  to  achieve  a 30°  diameter  circular  beam  are  supplied 
in  Figure  A.l.  (Actually,  the  amplitude  coefficients  that  were  achieved 
through  Goal  Programming  have  been  normalized  to  unity.  This  shading 
process  resulted  in  a gain  reduction  of  fifteen  decibels  at  boresight.) 

The  two-dimensional  contour  map  in  Figure  A. 2 represents  the  radiation  of 
the  transducer  array  as  a function  of  direction  (0,0).  The  response, 
measured  in  negative  decibels,  does  Indeed  provide  a 30°  diameter  circular 
beam  which  meets  or  surpasses  the  necessary  side  lobe  dropoff  throughout 
most  of  the  auxiliary  lobes.  Figure  A. 3 is  the  three-dimensional 
representation  of  the  radiation  of  the  transducer  array  with  these 
characteristics.  Although  the  perspective  of  this  plot  is  not  optimal, 
the  presence  of  the  circular  beam  is  still  apparent. 
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Figure  4.1.  Amplitude  and  Phase  Shading  Coefficients  for  a Circular  Beam 
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I Figure  4.2.  Contour  Map  of  the  Circular  Beam  Response  Pattern 
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4. A Beam  2 : A Rectangular  Array  Pattern 

Figure  A. A provides  the  amplitude  and  pltase  coefficients  which 
Goal  Programming  found  to  best  satisfy  the  72"  x 36"  rectangular 
specifications.  (Here,  again,  the  "optimal"  amplitude  coefficients 
were  normalized  to  unity.  In  this  case,  the  .shading  process  reduced 
the  gain  by  twenty-two  decibels  at  boreslght.)  The  contour  map  of  the 
response  pattern  thus  achieved  Is  found  In  Figure  A. 5.  The  only 
deviations  from  the  desired  response  levels  are  the  slight  violations 
found  at  the  corners  of  the  main  beam  and  the  side  lobe  levels  Immediately 
^adjoining  the  main  beam.  These  violations  are  minor  as  can  be  seen  from 
the  three-dimensional  plot  of  the  radiation  pattern  in  Figure  A. 6. 
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Figure  4.4,  Amplitude  and  Phase  Shading  Coefficients  for  the  Rectangular  Beam 
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Figure  4.5.  Contour  Map  of  the  Rectangular  Response  Pattern 
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CHAPTER  V 


I 
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SUMMARY  A-ND  CONCLUSIONS 
5 . 1 

The  broad  underlying  objective  of  this  study  was  to  demonstrate 
the  Goal  Programming  approach  in  the  analysis  and  solution  of  multi- 
criteria  system  design  problems. 
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The  particular  application  studied  has  been  the  design  of  planar 
arrays  in  order  to  achieve  specific  response  objectives.  The  formulation 
and  achievement  of  the  12  x 12  planar  array  designs  outlined  within  this 
thesis  provide  an  indication  of  the  typical  response  objectives 
encountered  and  the  ability  of  Goal  Programming  to  satisfy  the  multiple, 
conflicting  goals  of  such  designs. 

The  specific  objectives  of  this  thesis,  as  given  in  Chapter  I, 

were : 

1.  Develop  a mathematical  model  which  describes  the 
multicriteria  array  design  problem  via  a set  of 
mathematical  statements  and  extend  this  model  to 
the  goal  programming  formulation. 

2.  Achieve  two  specific  acoustic  array  designs  using 
the  Goal  Programming  approach. 

3.  Discuss  the  versatility  of  Goal  Programming  in 
array  design  relative  to  other  tools  used  in  such 
designs. 

4.  Through  the  successful  nature  of  the  results 
achieved  and  the  straightforward  nature  of  the 
problem  formulation,  encourage  the  use  of  Goal 
Programming  in  resolution  of  even  more  complex 
array  design  problems. 


The  first  objective  Is  responded  to  In  Chapter  III.  The  chapter 
on  empirical  results.  Chapter  IV,  provides  the  results  of  the  second 
objective.  Chapters  II  and  IV  combine  to  satisfy  the  third  objective. 
Finally,  it  is  hoped  that  the  achievement  of  the  fourth  objective  will 
be  encouraged  by  the  accomplishment  of  the  first  three. 

In  support  of  the  fourth  objective,  a summary  of  the  process 
necessary  to  use  the  computerized  Goal  Programming  approach  to  array 
design  follows: 

1.  The  development  of  a preliminary  mathematical  model 
which  essentially  describes  the  multicriteria  design 
problem  via  a set  of  mathematical  statements. 

2.  The  extension  of  the  preliminary  mathematical  model 
to  tlie  goal  programming  formulation.  This  involves 
the  addition  of  deviation  variables  in  the  objective 
functions  and  the  statement  of  the  "achievement 
function"  in  terms  of  the  deviation  variables. 

3.  The  development  of  the  subroutine  (YVALUE)  which 

defines  the  objectives  of  the  system  design  in  the 

FORTRAN  computing  language.  (A  listing  of  that 

» 

subroutine  for  the  design  of  the  circular  beam 
presented  herein  may  be  found  in  Appendix  C.) 

4.  The  solution  of  the  design  problem  using  subroutine 
YVALUE  in  conjunction  with  the  nonlinear  Goal 
Programming  code  as  listed  and  documented  in  the 
appendix  of  Reference  3.  The  use  of  this  code  is 
facilitated  by  the  user's  guide  presented  in 
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Reference  3,  as  well  as  the  analysis  of  the  procedure 
provided  in  Appendix  B of  this  thesis. 

5 . 2 Conclusions 

It  is  felt  that  the  basic  nature  of  most  real  design  problems 
should  strongly  emphasize  the  need  for  Goal  Programming  with  its 
capacity  to  resolve  multiple,  conflicting  design  objectives.  By  the 
nature  of  this  capacity,  the  Goal  Programming  model  has  no  constraint 
functions  and  thus  more  efficient  optimization  techniques  are  employed 
in  its  solution. 

Another  highly  favorable  characteristic  of  the  Goal  Programming 
technique  to  be  reiterated  in  conclusion  is  its  straightforward  extension 
to  other  complex  multicriteria  design  problems.  It  is  hoped  that  this 
easy  transition  to  more  complex  problems  coupled  with  the  success  of 
this  research  endeavor  will  stimulate  the  use  of  Goal  Programming  in 
further  system  designs. 


APPENDIX  A 


OUTPUT  FROM  AN  ARRAY  CONSISTING  OF 
SETS  OF  SYMMETRICALLY  LOCATED  ELEMENTS 


A.l  The  Derivation  of  the  Power  Response  Objectives  of  a Symmetric 
Planar  Array 

Consider  an  element  of  a planar  transducer  array  located  at  a 
point  whose  polar  coordinates  are  (r,6)  In  the  X-Y  plane  (Figure  4,  see 
text).  The  output  from  a signal  arriving  from  a direction  is: 


P . = Ae 

r,6 


j (o)t  + u) 


= Ae~^  (cos  u + j sin  u)  , 


(A.l) 


where ; 


u = kr  sin  0 cos(())  - 6)  + a 


= kx  sin  0 cos  4>  + ky  sin  9 sin  (J)  + a 


(A. 2) 


A is  the  shading  coefficient 
k is  the  wave-number  2it/X 


a is  the  element  phasing. 


The  output  of  an  array  of  M such  elements  is: 


’ * Z A e"^  (cos  u + j sin  u ) 

tot  , m m m 

m=l 


(A. 3) 


The  term  e'*  merely  indicates  the  frequency  dependence,  and  is  dropped 
when  the  equation  is  used  for  the  analysis  of  directional  response. 
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Taking  advantage  of  the  symmetry  inherent  in  the  problem,  the 
expression  for  the  total  output  can  be  further  simplified.  For  a pair 
of  elements  at  opposite  ends  of  a diameter  through  the  origin,  i.e.,  at 
(r,6)  and  (-r,5)  with  the  same  phase  shift  a applied,  we  have: 


P = + P = A[cos(u'  + a)  + cos(a  - u') 

+ j sln(u'  + a)  + sln(a  - u')] 
= 2 A cos  u'(cos  a + j sin  a)  , (A. 4) 


where  u'  = u - a . 

Thus,  the  output  of  an  array  of  M such  pairs  may  be  written; 


M 

P^  = 2 Z A cos  u' (cos  a + j sin  a ) 

btot  , m m m m 

m=l 


(A. 5) 


For  M sets  of  four  elements  In  quadrantal  symmetry  [whose 

cartesian  coordinates  may  be  written  (x  ,y  ) , (-x  ,y  ) , (x  ,-y  ) , 

mm  mm  mm 

(-X  ,-y  )]  all  with  the  same  coefficient  A and  phase  a , we  may 
m m m m 

incorporate  that  symmetry  to  obtain: 


P--  = 4 Z A (cos  01  + j sin  a ) cos(k  x sin  G cos  tp) 

ZQ  m m m m 

cos(k  y sin  G sin  ip) 
m 


If  we  are  not  Interested  in  the  phase,  we  may  write  for  the  power 
response; 


‘^ZQ  “ 


M 

Z A cos  a cos  U cos  V ) 
- tn  m in  m 

m=l 


M 

+ ( Z A sin  a cos  U cos  V ) 
• m m in  m 

m»l 


(A. 6) 


i 
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where ; 


APPENDIX  B 


THE  EXPERIMENT 

B. 1 Int  roduct Ion  to  the  Modified  Hooke  and  peeves  Algorithm 

In  order  to  describe  the  procedure  used  to  achieve  a "satisfactory" 
design.  It  Is  necessary  to  further  elaborate  on  Ignlzlo's  extension  of 
the  pattern  search  method  of  Hooke  and  Jeeves  employed  by  the  nonlinear 
Goal  Programming  algorithm  to  solve  the  multiple  objective  nonlinear 
model  [3]. 

Simply  stated,  pattern  search  begins  with  an  initial  base  point, 
provided  as  Input  to  the  algorithm,  and  perturbs  about  this  point  along 
all  the  coordinates.  In  this  manner,  it  searches  along  all  dimensions 
accelerating  the  search  If  previous  searches  have  been  successful  and 
decreasing  otherwise.  Input  to  the  algorithm  Includes  the  Initial  step 
sizes  for  each  dimension,  an  acceleration  factor  to  be  used  upon  success, 
and  a step  reduction  factor  to  be  employed  when  search  In  all  directions 
falls  to  Improve  the  solution.  Thus,  pattern  search  begins  at  the 
specified  Initial  point  and  searches  along  each  axis  at  the  step  size 
given.  If  it  succeeds  in  any  direction,  search  advances  from  this 
successful  point  along  the  next  coordinate.  After  searching  along  each 
dimension,  the  solution  vector  achieved  as  the  result  of  the  search  Is 
tested  to  see  If  It  gives  a preferred  achievement  vector.  If  the  new 
solution  has  an  improved  achievement  vector,  a , the  resultant  point 
becomes  the  new  base  point.  The  search  is  then  accelerated  by  renewing 
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the  search  at  a trial  point  found  at  the  end  of  a vector  from  the  old  base 
point  through  the  new  base  point  at  a length  of  the  acceleration  factor 
times  tfie  distance  from  the  old  and  new  base  points.  Therefore,  when  the 
search  succeeds  in  a given  direction,  pattern  search  assumes  that  the 
search  is  heading  in  the  right  direction  and  accelerates  along  this  path. 

If  search  about  a point  fails  to  improve  the  achievement  vector  values, 
pattern  search  returns  to  the  most  recent  base  point  (the  most  successful 
solution  to-date)  and  perturbs  around  it  at  reduced  distances  of  the 
reduction  factor  times  the  step  size.  Thus,  failure  indicates  that  the 
optimal  solution  does  not  lie  in  this  direction  or  has  been  overstepped 
and  pattern  search  attempts  to  backtrack  and  reduce  the  step  size  of  the 
. search.  Any  subsequent  failures  to  improve  the  solution  are  met  with 

further  reduction  of  the  step  sizes  by  further  multiplication  by  the 
reduction  factor. 

Termination  of  the  search  algorithm  will  occur  when  the  search  has 
achieved  a local  optimum  which  may  be  the  global  value.  In  order  to 
facilitate  recognition  of  this  optimum  and  avoid  the  continuation  of  the 
algorithm  when  it  is  only  improving  the  base  point  by  a very  small 
Increment  (often  a computational  round-off  factor),  the  algorithm  input 
specifies  a termination  factor.  Whenever  two  sequential  iterations  in  the 
search  lead  to  achievement  values  that  differ  by  no  more  than  this  factor, 
the  search  is  completed.  It  was  necessary  to  implement  another  rule  of 
termination  as,  oftentimes,  the  pattern  search  algorithm  can  continue 
advancing  the  solution  by  a small  design  Improvement  at  each  iteration. 

• This  improvement  when  weighed  against  the  computation  time  involved  in 

each  iteration  may  not  be  sufficient  to  warrant  continuation  of  the  search. 

f 
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The  rule  provided  is  to  terminate  the  search  algorithm  after  a prespecified 
number  of  Iterations.  This  is  anotlier  factor  which  is  provided  by  the 
user  for  the  Goal  Programming  code.  Another  safeguard  provides  against 
pattern  cycling  upon  successive  failures  and  is  a user  designated  factor 
which  dictates  the  maximum  number  of  Increment  reductions. 

With  this  brief  introduction  to  the  pattern  search  algorithm, 
several  factors  have  been  brought  to  light.  Input  to  the  algorithm  must 
provide  the  initial  starting  point,  the  initial  search  step  sizes,  the 
acceleration  factor,  the  step  reduction  factor,  the  maximum  number  of 
search  patterns,  and  the  maximum  number  of  Increment  reductions. 


B.2  Algorithm  Starting  Points 


To  compensate  for  the  complex,  multimodal  nature  of  antenna  pattern 
response  functions,  each  design  problem  formulation  has  been  subjected  to 
several  solution  attempts  using  dlftcrent  starting  points.  It  is  apparent 
that  the  closer  the  starting  point  is  to  the  optimal,  the  more  quickl> 


the  algorithm  is  likely 
iteratively  repeating  a 
of  parametric  changes  to 
design  arrived  at  in  the 
for  continued  analysis, 
problem  formulations,  it 
points  be  evaluated  with 


to  converge.  Keeping  this  point  in  mind,  while 
specific  array  design  problem  with  any  of  a number 
Improve  upon  the  pattern  acliieved,  the  antenna 
prior  Iteration  was  used  as  the  initial  point 
In  an  attempt  to  save  computation  time  in  larger 
is  suggested  that  many  randomly  generated  starting 
the  "best"  of  these  employed  as  input  to  the 


algorithm. 


kp. 


i 
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B.3  Al_gorlthm_  Init  lal  Step  Sizes 

The  Initial  decision  on  step  sizes  is  an  intuitive  one.  If  the 
optimal  solution  Is  thought  to  be  relatively  close  to  the  specified 
starting  point,  the  search  should  begin  at  small  Intervals.  The  choice 
of  the  "best"  step  sizes  is  at  first  an  Iterative  procedure,  determining 
those  step  sizes  which  consistently  overstep  or  underachieve  the 
Increment  necessary  to  achieve  the  optimal.  Use  of  the  appropriate  step 
sizes  minimizes  the  computation  time  necessary  by  minimizing  the  number 
of  step  reductions  or  accelerations  necessary.  To  elaborate,  if  the 
Initial  step  size  is  too  large,  the  algorithm  may  have  to  undergo  a 
series  of  reduction  cycles  before  successfully  moving  to  an  improved 
point  in  a given  dimension. 

B . 4 Algorithm  Acceleration  and  Reduction  .ctors 

The  choice  of  acceleration  and  reduction  factors  is  also  based  on 
exposure  to  a number  of  different  values.  The  optimal  value  for  both  are 
characteristic  of  the  problem  at  hand  as  the  factors  are  dependent  upon 
the  functions  to  be  optimized.  Any  inadequacies  in  the  acceleration  or 
reduction  factors  are  compensated  for  by  the  iterative  process  of 
acceleration  and  reduction,  but  again,  the  computation  time  consumed  by 
this  process  can  be  alleviated  by  adjusting  these  values. 

B . 5 Achievement  Function  Tolei'ances 

The  achievement  function  tolerances  are  also  dependent  upon  the 
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characteristic  values  of  the  achievement  function  vector.  In  some  design 
problems,  a difference  of  0.01  in  achievement  function  values  is 
significant,  whereas,  in  others  this  is  not  an  accurate  or  valuable 
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measure  of  Improvement.  The  larger  the  values  of  the  tolerances,  the 
tighter  the  default  constraint  on  achievement  function  Improvements; 
thus,  rapid  termination  of  the  algorithm  can  be  brought  about  by  setting 
high  tolerances.  This  measure  is  of  use  when  computation  time  becomes  a 
serious  consideration  due  to  the  size  of  the  problem.  By  means  of  the 
specification  of  the  maximum  number  of  search  patterns  and  the  maximum 
number  of  Increment  reductions,  two  additional  options  are  available  to 
regulate  the  computation  time. 


B.6  The  Number  of  Response  Objectives 

The  initial  decision  concerning  the  number  of  response  functions 
to  include  as  the  design  objectives  (functions  of  0 and  ) is  also 
basically  an  intuitive  one.  Initially,  the  number  of  response  objectives 
should  be  determined  on  the  basis  of  any  previous  knowledge  or  intuition 
concerning  the  ease  of  attaining  the  desired  pattern.  If  the  pattern  is 
naturally  "easy"  to  attain,  the  number  of  response  functions  defined  in 
the  objectives  may  be  few.  If  the  desired  response  is  an  imposing  one, 
numerous  objectives  may  have  to  be  Included  in  the  model  in  order  to 
control  response  inadequacies  that  would  occur  within  sizeable  uncon- 
strained regions.  It  is  best  to  begin  with  a conservative  estimate  of 
the  number  of  objectives  necessary  as  this  allows  the  potential  to  see 
where,  if  at  all,  the  response  pattern  is  naturally  weak. 

In  choosing  the  0,i))  values  to  be  Included  in  the  response  objec- 
tives, it  is  best  to  begin  with  equally  spaced  radiation  levels  as  unequal 
spacing  may  apply  an  unnecessary  amount  of  pressure  on  the  more  closely 
constrained  area  of  response  at  the  expense  of  the  radiation  levels  In 
the  region  with  few  constraints. 
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If  upon  inspection  of  the  pattern  achieved  by  this  Initial  attempt 
it  is  found  that  an  unconstrained  area  of  the  response  fails  to  meet  the 
desired  level  and  is  possibly  even  affecting  tlie  response  within  the 
constrained  region,  the  "grid"  of  the  objectives  In  this  area  should  be 
increased.  This  step  should  be  taken  with  tlie  understanding  that  the 
radiation  level  or  response  may  merely  shift,  although  the  effects  may  be 
diffused  into  a larger  area  causing  a problem  not  nearly  as  drastic  as 
the  one  it  alleviated.  If  many  undesirable  response  levels  occur  within 
the  unconstrained  regions,  the  logical  step  is  to  tighten  the  grid  of 
design  objectives.  In  other  words,  an  ovei.ili  weak  pattern  due  to  tlie 
response  in  unconstrained  response  regions  may  generally  be  alleviated  by 
Increasing  the  radiation  levels  (functions  of  6,<}i  ) Included  in  the 
objectives  of  the  model. 

B . 7 We ijilit  i n^  F acto r s 

Unless  some  Insight  is  available  concerning  the  response  difficulties 
of  a particular  design  before  building  the  model,  it  Is  best  initially  to 
give  equal  weight  to  the  design  objectives  of  both  the  main  beam  and  side 
lobes.  The  pattern  arrived  at  when  the  nonlinear  Coal  Programming 
algorithm  attempts  to  fulfill  these  equally  weighted  objectives  will  give 
some  indication  of  the  Inherent  difficulties  in  the  design  of  the  desired 
pattern.  If  the  pattern  achieves  or  exceeds  the  design  objectives  on  the 
side  lobes  while  the  main  beam  or  regions  within  the  beam  fall  to  meet  its 
desired  levels,  the  model  should  be  restated  with  additional  weights  on  the 
main  beam  objectives,  forcing  the  algorithm  to  place  additional  importance 
on  the  response  levels  at  the  designated  areas  of  the  beam.  The  decision 
as  to  the  magnitude  of  these  weights  is  gained  through  actual  experience. 


51 


b 


¥ ‘ 


t 


I 


4 


with  numerical  values  determined  by  the  degree  of  difficulty  characteristic 
of  the  region  to  be  weighted. 

As  was  previously  stated,  the  arrival  of  the  "correct"  weighting 
factors  for  a given  design  is  often  an  iterative  process  with  the 
completion  of  the  process  signaled  upon  arrival  at  a "satisfactory" 
antenna  array  design. 

B . 8 The  Response  Levels 

When  establishing  the  levels  to  be  achieved  by  the  design  objectives, 
the  difficulty  or  Impossibility  of  achieving  a particular  pattern  at  the 
levels  specified  may  not  be  known.  It  may  be  found  that  unrealistic 
goals  have  been  set,  e.g.,  the  side  lobes  cannot  maintain  an  overall 
response  level  as  far  removed  from  that  of  the  main  beam  as  specified. 

A judgment  of  this  nature  should  not  be  made  without  evaluating  the  results 
of  a number  of  starting  points,  a range  of  weighting  factor  configurations, 
and  a sufficient  number  of  objectives.  If  it  Is  found  that  unrealistic 
goals  have  been  set.  It  Is  possible  to  Iteratively  reformulate  the  model 
and  establish  an  upper  bound  on  the  response  levels  that  can  be  achieved 
for  a particular  response  pattern. 

The  process  of  adjustment  can  begin  by  bringing  down  the  specifica- 
tions for  an  entire  region  to  levels  easily  achieved  by  most  of  the  region. 
The  algorithm  should  then  focus  Its  attention  on  the  achievement  of  these 
adjusted  levels  In  the  trouble-spots  as  all  other  radiation  levels  have 
achieved  the  goals  set.  If  the  pattern  responds,  then  the  specifications 
can  be  iteratively  Increased  seeking  the  upper  limit  on  the  response  that 
can  be  achieved. 

This  concludes  the  analysis  of  the  procedure  used  in  the  design  of 
planar  antenna  arrays  with  the  computerized  nonlinear  Goal  Programming 


algorithm. 


APPENDIX  C 


SUBROUTINE  YVALUE  TO  PROVIDE  THE 
MATHEMATICAL  MODEL  FOR  THE  12  x 12  ARRAY 
ACHIEVING  CIRCULAR  BEAM  SPECIFICATIONS 
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